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Abstract 

We investigate the physical interpretation of the Riemann zeta function as 
a FZZT brane partition function associated with a matrix/gravity correspon- 
dence. The Hilbert-Polya operator in this interpretation is the master matrix 
of the large N matrix model. Using a related function 2(z) we develop an 
analog between this function and the Airy function Ai(z) of the Gaussian ma- 
trix model. The analogy gives an intuitive physical reason why the zeros lie 
on a critical line. Using a Fourier transform of the S(z) function we identify a 
Kontsevich integrand. Generalizing this integrand tonxn matrices we develop 
a Kontsevich matrix model which describes n FZZT branes. The Kontsevich 
model associated with the 5(z) function is given by a superposition of Liouville 
type matrix models that have been used to describe matrix model instantons. 

1 Introduction 

It is an old idea that if a Hermitian operator can be found which has eigenvalues of 
the form A n = —i(p n — \) where p n are the nontrivial zeros of the zeta function then 
the Riemann hypothesis would be true. This would follow because the eigenvalues 
of a Hermitian operator are real. No such operator has yet been found however. 

In the Heisenberg matrix formulation of quantum mechanics [1] one represents 
observables with infinite matrices which are Hermitian. The eigenvalues of the 
matrix are what are measured in an experiment and hence are real. The infinite 
matrix can be constructed by forming an N x N matrix and taking the large N 
limit. No such large N matrix whose eigenvalues are related to the Riemann zeros 
has been found. 

One can also consider theories called matrix models in which the dynamical 
variables are such large N matrices [2] [3] [4]. There is a remarkable correspondence 
between such matrix theories and continuous theories describing a quantum theory 
of world sheet gravity and low dimensional string theory. In such a correspondence 
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invariants of the matrix theory are related to geometric observables in the world- 
sheet gravity usually through an integral transform [4]. This integral transform 
takes one from a variable in the expansion of the characteristic polynomial of the 
large N Hermitian matrix description to a Liouville variable describing the size of 
the string or of a of a macroscopic loop of 2d gravity in the continuum description. 
The correspondence arises because Feynman graphs in the matrix theory description 
can yield discrete representation of surfaces which become continuous as one takes 
N to infinity [5]. 

In this paper we interpret the Riemann zeta function as being related to a 
particular observable in the matrix/gravity correspondence namely the FZZT brane 
partition function of a matrix model and interpret it's master matrix as the Riemann 
operator. The potential for this matrix model is more complicated than most of the 
ones considered in the literature. Nevertheless the techniques of simpler matrix 
theories can be applied to this case as well. This paper is organized as follows. 
In section 2 we discuss the master matrix approach to matrix models. We discuss 
some of the conceptual advantages of the approach as well as the difficulties. In 
section 3 we discuss introduce the FZZT brane partition function from the matrix 
model point of view. In section 4 we determine the Kontsevich integrand associated 
with the Riemann zeta function and develop an analogy between the Riemann zeta 
function and the Airy function which is the FZZT partition function of the (2, 1) 
minimal matrix model. In section 5 we discuss how to approximate the matrix 
model associated with Riemann zeta function using the generalized (p, 1) matrix 
model for large p whose FZZT partition function is a generalized Airy integral. In 
section 6 we review the main conclusions of the paper. 



2 Master matrix 

If one can find a special infinite Hermitian matrix Mq such that: 

E(z) = det(M - zl) 

where 

m = c(zz + \)v^ + l )7r -vv-*/ 2 (-^ _ I) 

then the Riemann hypothesis would be true. This is because this function can be 
written in product form as: 
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The eigenvalues of the Hermitian matrix Mq are denoted by \ n and are related to 
the Riemann zeros via p n = i\ n + 1/2. Then the product becomes: 

"( \ - Irrn _ iz + 1 / 2 \ _ 1 tt A » ~ z 

a W- 2 ll 11 iA n + l/2 , "2 i n i A n -»/2 

This vanishes at the values A ra just as the formal determinant expression. The \ n 
are real if the matrix Mq is Hermitian and thus the Riemann Hypothesis would be 
true. Unfortunately just as for the Riemann operator referred to above no such 
infinite matrix Mq has ever been constructed. 

The difficulty in constructing Mq is somewhat similar to the difficulty in con- 
structing a master field or master matrix in large N field theory of matrix theory 
[6] [7]. A master field or master matrix is a special large N matrix such that statisti- 
cal averages of an observable can be computed by simply evaluating the observable 
on the the special large N matrix. The reason that a master matrix exists is because 
at large N expectation values factorize as: 

= (0!)(0 2 ) + 0(l/iV 2 ) 
((0-(0)) 2 ) = (0 2 )-(0) 2 = 0(l/iV 2 ) 

where: 

(O) = J DMO(M)e- v ^ 

and V(M) is a matrix potential. Thus variances vanish so the observable's value is 
localized on a particular matrix as N — > oo just as particle trajectories are localized 
on classical solutions as h goes to zero. Once such a master field is found the above 
observables are simply given by: 

(O) = O(Mq) 

There are several such observables in matrix theory. We discuss some of these in 
the next section. 

For a general matrix model with potential V(M) the master matrix can be 
written [6] [7]: 

oo 

M = S- l TS = S-\a + t n a +n )S 

n=0 

where the similarity transformation S is defined so that Mq is Hermitian and the 
operators a, a + obey [a, a + ] = I. One can expand the master matrix as a function 
of the Hermitian operator x = a + a + as: 

Mq{x) = gix + g 2 x 2 + ... 
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One can also define an associated complex function: 

-i oo 

M (y) = - + E f ny n 

as well as a conjugate matrix Po that satisfies: 

[Po,M ]=I 

The Master matrix can be determined from the equation [6] [7] : 

(V(M (x )) + 2P ) |0) = 

Here |0) is the vacuum state annihilated by a. The master matrix is closely connected 
with the resolvent R{z) and eigenvalue density p{x) through: 

C 

The associated function Mo(y) obeys the relation: 

R(M (y)) = M (R(y)) = y 

The function yMo(y) is the generating functional of connected Green functions for 
the generalized matrix model. While the concept of the master matrix is appealing, 
to construct the master matrix explicitly is equivalent to finding all the connected 
Green functions which amounts to solving the theory. This can be done for the 
potential V(M) = Tr(M 2 ) but for the general matrix model is quite difficult. In the 
next three sections we turn to other methods of dealing with the generalized matrix 
model which are somewhat more tractable and apply them to the interpretation of 
the zeta function. 

3 FZZT brane 

One observable of matrix models is the exponentiated macroscopic loop or FZZT 
brane partition function [8] [9] [10] [11]. This is given by: 

B(z) = det(M - zl) 

This is the characteristic polynomial associated with the matrix M. It's argument 
z can be complex. In the context of the Riemann zeta function Q{s) the variable is 
related to the usual argument of the zeta function by s = iz + \. Another observable 
is the macroscopic loop which is the transform of the Wheeler-DeWitt wave function 
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defined on the gravity side of the correspondence [12]. 

oo 

W(z) = -Tr\og{M - zl) = lim( / ^Tr{e^~ zI+M ^>) + loge) 

£ 

where e is a UV cutoff. 

The resolvent observable mentioned above is defined by: 

Finally one has the inverse determinant observable defined in [12]. 

If a special master matrix Mq can be found then expectation values such as 

(B(z)) = (det(M - zl)) = J DMdet(M - zI)e' v{M) = det(M - zl) 

reduce to evaluating the observable at Mq. In the context of the E(z) function the 
desired relation is of the form: 

E(z) = det(M - zl) = (B(z)} = (det(M - zl)} = J DMdet(M - zI)e' v{M) 

Some matrix potentials that have been considered are 

V(M) = Tr(M 2 ) 

which describes 2d topological gravity or the (2,1) minimal string theory [13] [14] 
[15]. A quartic potential: 

V(M) = Tr{—M 2 + gM A ) 

is used to describe minimal superstring string theory [16] [17] [18] [19]. A more com- 
plicated matrix potential is 

V(M) = —Tr(M + Iog(J - M)) = V —Tr(M m ) 

m=2 

which defines the Penner matrix model [20] [21] [22] [23] [24] and is used to compute 
the Euler characteristic of the moduli space of Riemann surfaces. 

Another matrix model that has been introduced is the Liouville matrix model 
[25] [26] with potential given by: 



V(M) = Tr(aM + fie 
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with cosmological constant \i so that: 

e -V(M) _ e -oTrM e -iiTre M 

In this paper we will encounter the matrix potential determined by: 

OO Q 

e -U(M) = g ((? V e 2TrM _ t^^TrMy-^Tr^) ^ 
9=1 2 

The partition function for this matrix model can be seen as a superposition of 
partition functions of Liouville matrix models with cosmological constants of the 
form. 

H = q 2 ir 

for integer q. The origin of this particular matrix model and it's relation to the zeta 
function will be discussed in the next section. 

4 Kontsevich integrand 

To see how the matrix potential (3.1) arises it is helpful to consider how the co- 
efficients of the characteristic polynomial observable B{z) can be determined by 
expanding as a series in z. If the function is interpreted as a characteristic 

polynomial then one can obtain these coefficients from the expansion: 

oo 

i 

OO o 

9=1 2 

Inserting the coefficients a,2 n into and summing over n we can represent 
as an integral transform: 

oo oo 

E(z) = 4 J ^(«+V2)/2 J2 (q\ 2 £ 2 - ^q 2 7r£)e-^ = 4 J ^(«+i/2)/2^i/2 



where 
and 
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Defining the variable (/> by £ = we have: 

r. CO Q 

E[z] = / £ (ttW 1 * - ^V^"** 2 '* (4.1) 

fc=i 

which is a well known integral expression for the function 

For the simple potential V(M) = Tr(M 2 ) the exponentiated macroscopic loop 
observable (FZZT brane) can be computed. It is given by the Airy function [15]: 

Ai(z) = J DMdet(M - zI)e~ Tr ^ = J #e i2 ^H (4.2) 

Because this function is associated with an Hermitian matrix model it's zeros are 
real. This is the analog of the Riemann hypothesis for V(M) = Tr(M 2 ). The 
similarity between the integral representations of (4.1) and (4.2) suggest an analogy 
between the Airy and zeta functions. 

To illustrate a comparison between the Airy function and the zeta function 
consider figures 1 and 2. The zeros disappear as one moves off the critical line which 
corresponds to z real in both cases. This suggests that the zeta function corresponds 
to a Hermitian matrix model. Table 1 illustrates the comparison on both sides of the 
correspondence. The question mark indicates the (substantial) missing information 
involved in a matrix/gravity approach to the Riemann hypothesis. 

Qualitative differences exist between between the functions Ai(z) and The 
Ai(z) function is exponentially decaying to the positive z axis. This is a result 
of Stokes phenomena where an exponentially growing form of the Airy function 
is completely absent in the right z axis. For the H(z) function one does not see 
exponentially decaying function in the positive z axis. Instead one has identical 
behavior in the positive and negative z axis. One way to see the difference is to use 
a Riemann-Hilbert Problem approach to both functions [27] [28] [29] [30] [31] [32]. In 
the case of the Airy function this leads to the differential equation [27]: 

Ai"(z) = zAi(z) 

whereas in the case of the H(z) function one does not obtain a differential equation 
but a discrete equation [27] : 

E(z) = E(-z) 

Indeed it is known that the zeta function does not obey a finite order differential 
equation so this may be a possible explanation for the qualitative difference between 
the two functions. It would be interesting to explore further the differences between 
the two functions using the Riemann-Hilbert approach of and their interpretations 
as FZZT brane partition functions. . 

The integral representation of the Airy function has a matrix integral general- 
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Figure 1: Magnitude of the Airy function on the real axis. The zeros are all located 
on the negative real axis. An intuitive way to understand this is that the Airy 
function is the FZZT brane partition function of a matrix model with potential 
V = Tr(M 2 ) and Kontsevich integrand e~ u ^ = exp(i^> 3 /3) 




Figure 2: Magnitude of the function Q(ix + 1/2) on the real axis. The zeros are 
all symmetrically located on the real axis. An intuitive way to understand this is 
that the S function is the FZZT brane partition function of a matrix model with a 
suitably chosen potential V(M) and Kontsevich integrand e~ u ^ . 
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Observable 


General 


Airy 


Zjeta 


Master Matrix 

Potential 

FZZT Brane 
Macroscopic Loop 

Kontsevich Integrand 


M 

V(M) 

B{z) 
W(z) 

e -U(4>) 


a + a + 

Tr(M 2 ) 

Ai(z) 
log Ai(z) 


? 

lim Tr{V p (M) + P f: s k V k (M)) 

p— >oo ^ ^ 

C(*2 + l/2) 
IogC(w+l/2) 

£ (gVe 2 * - §g 2 7re^)exp(-7rgV) 

9=1 



Table 1: Analogy between the Airy function and the Riemann zeta function. The 
quantities Vjfe(M) and defined by a generalized (p, 1) matrix model in the following 
section. 

ization. The matrix potential is defined from: 

e -U(9) = e i|Tr(* 3 ) 

The matrix generalized Airy function is given by: 

Ai(Z) = J d*<F r We- u W 

In the above and Z are n x n matrices. The interpretation of this matrix integral 
is that it describes n FZZT branes. The matrix $ in the Kontsevich integrand is 
an effective degree of freedom describing open strings stretched between n FZZT 
branes [33] [34] [35] [36] . 

One can try to interpret the integrand of the function in a similar manner. 
In that case the analog of the potential defined by: 

OO r, 

e -U(*) = (7r 2 fc 4 e 2Tr« _ f^gTrijg-^Tre* 
fc=l 2 

and the analog of the the matrix integral describing n FZZT branes is: 

/OO o 
D$ e «V(Z*) g (7r 2 fc 4 e 2Tr« _ f^gTrS^-^Tre* 

fc=l 2 

This is the origin of the matrix model given by (3.1). As discussed in section 2 this 
can treated as a sum of Liouville type matrix models. 

5 Relation to generalized (p, 1) matrix models 

The Airy function is the FZZT partition function for the (2, 1) minimal matrix 
model. In [37] the FZZT partition function was given for the generalized (p, 1) min- 
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imal matrix model with parameters s&. This theory has a characteristic polynomial 
or FZZT partition function given by: 



B(z) 



1 

2^ 



be i Z <t>- ^ (i4>r +1 +El=! s * fft w) fc+1 



Unlike the (2, 1) matrix model the definition of the generalized (p, 1) matrix model 
requires a two matrix integral of the form [37] [38] [39] [40] [41] [42] [43]: 



Z(p,i)(9) = J DMDAe~ s ( V ( M + I )~ AM ) 



Comparison with the integral representation of the E(z) function shows that a gener- 
alized matrix model for large p can be constructed as an approximation. This can be 
compared with the formulas from the previous section to compute the corresponding 
coefficients s k - One writes: 



log ( ]T (ir 2 ^ 2 * - -7rfcV)e 



-TTfcV 



P + 



fc+i 



In the above formula the function on the left is expanded to order p+1 in the variable 
(j). We denote this terminated expansion by E p (z). Another way to compute the 
coefficients s^is to differentiate the left hand side and set: 



— (fc+l) / oo o 

r = —r^~ d\ log £ W ^ ~ 
\fc=i 1 



■irk 2 e<t> 



From the integral representation one has: 



where: 



QEp(z) = zEp(z) 
PE p (z) = -d z E p (z) 



p-i 

k=0 



Inserting this operator into the above equation one has the generalization of the 
Airy equation given by: 



(PV + Y^s k P k )Z p (z)=zE p (z) 

k=0 



(5.1) 



To recover the equation for the full E(z) function one has to take p to infinity which 
agrees with the fact that the zeta function does not obey a finite order differential 
equation. 

Note that z and <p are in some sense canonically conjugate [37]. Denote the 
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Fourier transform of the E(z) function as S(p) then: 

The generalized Airy equation them becomes in Fourier space: 

p-i 

{<p + Y,*k<l> k )%>{<i>) = Q%>{<l>) 

k=0 

This can be written: 

(U'((f>) - Q)E(<t>) = (5.2) 

where: 

OO Q 

fc=i 2 

Equation (5.2) is very similar to the equation for the master matrix (2.1). Indeed if 
we set: 

<P = M (y) 

z = P (y) 

we see that y can be thought of as coordinates of a parametrization of the Riemann 
surface M V) \ which is determined from the <j) and z constraint U'((f>) — z = 0. If we 
make these variables into operators through: 

(t> = M {a,a + ) 
z = P (a,a+) 

this classical surface is turned into a quantum Riemann surface similar to those 
studied using noncommunative geometry [44]. 

Once one has obtained the coefficients Sk one can define matrix potential asso- 
ciated with a finite N theory as [37]: 

p-2 

V(M) = \\mTr{V p {M) + £ s k V k (M)) 
p ^°° k=i 

where: 

V k (M) = f2~(M j -I) 

3=1 3 

This is the matrix potential of Table 1 in the previous section. 

A set of orthogonal polynomials with this matrix potential through the integral 
equation: 

p-2 

n \ r - lim (V p (y+1)+J2 s k V k (y+l))+2zy 1 
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Or equivalently though the generating function definition: 



p-2 

-E- 

* =1 - > ff-m 



lim (^(^+1)+^; Sk V k (y+l))+2zy °° 

n=0 



nl 



These are the generalizations of the integral and generating function definitions of 
the Hermite polynomials associated with the (2, 1) minimal model. 

We note that some other entire functions can be treated in a similar manner. For 
example the reciprocal factorial function j^jy = r(z+i) has a product representation: 



1 oo 

e" TT ll + -)e-' /n 



no n ( i+„ 

n=L 



and integral representation: 



U(z) 



be iz4 >e- e 



The identity 



implies the equation 



or: 



U(z - 1) II(z) 
1 1 



e~ 9 



U(z) U(z) 



This is the analog of equations (5.1) for the function. The product repre- 
sentation shows that the zeros of the inverse factorial function are of the form 
X n = —1, —2, —3, — The inverse factorial function is similar to the zeta function 
in that it does not obey a finite order differential equation. It is similar to the Airy 
function in that it has all it's zeros on the negative real axis. It differs from both 
the Airy and zeta function in that it's zeros are of a simple form namely the nega- 
tive integers. The integral representation of the reciprocal factorial function seems 
related to the Liouville matrix model with Kontsevich integrand e~ u ^ = e~ e<l> and 
generalized (p, 1) matrix model with Sk = h- The matrix integral representation of 
the Gamma function in terms of the Liouville matrix model has been discussed in 
[25]. 

Most of our analysis has centered on the matrix side of the matrix/gravity corre- 
spondence. The gravity side is related through an integral transform. For example 
the macroscopic loop observable associated with the Riemann zeta function is given 
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by: 



oc 



logC(»« + l/2) = / r iz - 1/2 W{£)dl 



o 

In terms of the \ n this observable takes the form [45]: 



W{i) = 



\og£ 



1 



£ 



n 



2cos(A n log I) 
W 2 \og£ 



£(£ 2 - l)log£ 



1 



The indefinite integral of this Wheeler-DeWitt wave function is connected to the 
prime numbers p through: 



The FZZT brane partition function can also be represented by prime numbers as: 



Both of the above formulas follow from the Euler product formula of the zeta func- 
tion. Much of the physical intuition about the meaning of the FZZT brane and the 
Wheeler-DeWitt wave function occurs on the gravity side of the correspondence. 
Thus the connection of number theory and gravity in this context is quite intrigu- 
ing. 

Finally to approach the generalized Riemann hypothesis using the 
matrix/gravity correspondence one can replace the Kontsevich integrand e~ u ^ 
with a modular function. Indeed such modular functions already arise in the ma- 
trix/gravity CFT2/ACIS3 correspondence between two dimensional conformal field 
theory and three dimensional gravity with negative cosmological constant [46]. 

6 Conclusion 

In this paper we have examined the Riemann zeta function as a FZZT brane partition 
function involved in matrix models. The FZZT description gives rise to the physical 
interpretation of the Riemann hypothesis, that the E(z) is an entire function and has 
zeros on the critical line with z on the real axis (this corresponds to the Re(s) = i). 
The zeros are interpreted as eigenvalues of the master matrix. The macroscopic loop 
observable and resolvent also have physical interpretations in terms of the matrix 
model. In the gauge gravity correspondence the macroscopic loop is identified with 
the Wheeler-DeWitt wave function of the 2d world sheet gravity. The variable z is 
identified with the boundary cosmological constant in the 2d gravity. The matrix 



I 




logC(**+l/2) = ££-p 



n(iz+l/2) 



p n 
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gravity correspondence is the mapping between the matrix quantities and the 2d 
gravity computations. In a string theory context these in turn describe target space 
time processes. 

The Kontsevich integrand was identified using the Fourier transform of the 
function. Replacing the z variable by n x n matrix Z and the Kontsevich integrand 
by a matrix integrand one obtains representation of a matrix model describing n 
FZZT branes. The Kontsevich integrand is given by a superposition Liouville matrix 
models that have been used to represent instanton matrix models for the c = 1 string. 

Some long standing issues are indicated in Table 1. To identify and interpret 
the Master matrix associated with the Riemann zeta function. 
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